Decoupling of self-diffusion and structural relaxation during a fragile-to-strong 
crossover in a kinetically constrained lattice gas 
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We present an interpolated kinetically constrained lattice gas model which exhibits a transition 
from fragile to strong supercooled liquid behavior. We find non-monotonic decoupling that is due 
to this crossover and is seen in experiment. 



The viscosity of glassformers increases by several or- 
ders of magnitude upon cooling near the glass transition 
temperature. For strong glassformers, this increase is 
Arrhenius, that is, the viscosity increases exponentially 
with inverse temperature, whereas, for fragileglassform- 
ers, this increase is super- Arrhenius 0, 0> U Recent 
experimental and theoretical evidence indicate that frag- 
ile glasses can undergo a crossover to strong behavior at 
low enough temperatures 0,0,0- Here, we show that 
a fragile-to-strong crossover can be modeled by a union 
of two kinetically constrained lattice gas models 0, , 
and that, for this model, the variation of the product of 
the self-diffusion constant and the structural relaxation 
time is non-monotonic with temperature as a result of 
the crossover. The viscosity of salol seems to exhibit 
a fragile-to-strong crossover with lowering temperature 
0, llfjj l , and non-monotonic decoupling behavior similar 
to what we find in our model has been observed for the 
product of the shear viscosity and dielectric relaxation 
time in salol [Tl|. 

Our interpolated model is based upon the kinetically 
constrained triangular lattice gas (TLG) models intro- 
duced by Jackie and Kronig y|. These two-dimensional 
models are variants of lattice models proposed by Kob 
and Andersen . Each site of the triangular lattice has 
six nearest neighbor sites and can hold at most one par- 
ticle. A particle at site i is allowed to move to a near- 
est neighbor site, i', if (i) i' is not occupied and (ii) the 
two mutual nearest neighbor sites of i and i' are also 
empty. These rules coincide with a physical interpreta- 
tion of steric constraints on the movement of hard core 
particles in a dense fluid @ . We call the model with these 
rules the (2)-TLG because both mutual nearest neigh- 
bors need to be empty in order to facilitate movement. 
The model where the constraints are more relaxed we 
call the (l)-TLG: movement is allowed as long as either 
of the mutual nearest neighbors is empty. As with other 
kinetically constrained lattice gas models, the TLG has 
no static interactions between particles other than those 
that prohibit multiple occupancy of a single lattice site. 
These models exhibit all the salient features of glassform- 
ers (see, for example, Q, 0, 0, such as dynamical 
slowdown, heterogeneity and lengthscale growth. More- 
over, the (2)-TLG is a fragi le g lassformer and the (1)- 
TLG is a strong glassformer [T^] ■ 

To study crossover behavior, we introduce a model 
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FIG. 1: Simulated structural relaxation time, r a (top), and 
self-diffusion constant, Ds (bottom), as a function of inverse 
temperature at various values of the asymmetry parameter, b: 
6 = (circles), b = 10~ 2 (stars), b = 10" 3 (pluses), b = 10~ 2 
(triangles) and 6=1 (squares). The same labelling conven- 
tion is used for all figures. Dashed lines are guides to the eye 
and, unless otherwise indicated, statistical uncertainties are 
the size of the symbols used. 



which interpolates between the (2)-TLG and the (1)- 
TLG. A similar model which interpolates between 
the symmetrically and asymmetrically constrained Ising 
chains (i.e., the Fredrickson- Andersen 0] an d East 
models) predicts crossover behavior of relaxation times 
p|. For our studies, we introduce a parameter b, rang- 
ing from zero to one, which controls the dynamics in the 
following way. When b = 0, the only allowed move set 
is that of the (2)-TLG. When b > 0, a move from the 
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(l)-TLG move set is allowed with probability b. Finally, 
when 6=1, the model is the (l)-TLG. Results for the 
extreme cases of b = and 6=1 have been presented 
elsewhere .12]. In the present study, we perform addi- 
tional simulations at three intermediate values of b and 
over several densities, p. At each (p, b) state point, 32 
to 64 independent trajectories of length 10-100 times the 
structural relaxation time, r Q , were collected and ana- 
lyzed. For these simulations, t q is the time for the self- 
intermediate scattering function at wavevector tt to reach 
1/e of its initial value. We define a reduced temperature, 
T, as — ln(l — p) = 1/T. For additional computational 
details, see [T^ |. 

The structural relaxation time in the fragile (2)-TLG 
scales as a double exponential in inverse temperature and 
the strong (l)-TLG scales as a single exponential in in- 
verse temperature 0, . We therefore anticipate the 
following behavior for the crossover model. For a value 
of b between and 1, we expect the system to relax as 
if it were the fragile TLG at high temperatures. This 
is because the time needed for relaxation via the fragile 
TLG dynamics is faster than waiting a time proportional 
to 1/6 to relax via strong dynamics, despite the fact that 
strong relaxation is faster. At some point as temperature 
is lowered, the time required to wait for a strong move be- 
comes comparable to the fragile relaxation time and the 
system crosses over to strong behavaior. This behavior is 
clearly seen in FIG.^ F° r every value of b between and 
1, the structural relaxation time initially follows the (2)- 
TLG curve at high temperatures and then, at some lower 
temperature, T x , it crosses over to (l)-TLG behavior. As 
we would expect, T x increases, as 6 increases. 

We now turn to relaxation at larger lengthscales where 
a similar crossover behavior is observed in the self- 
diffusion constant, D$- See FIG. ^ In a dynami- 
cally homogeneous system, one would expect the prod- 
uct D s r a to be a constant independent of temperature. 
That is, one expects mean field relations such as the 
Stokes-Einstein (SE) relation to be valid. In dynamically 
heterogeneous systems such as glasses, these mean field 
relations are dramatically violated: relaxation at long 
lengthscales decouples from relaxation at shorter length- 
scales 0,0,0,0. In fragile glassformers, D s r a has 
been found to increase several orders of magnitude as 
temperature is lowered below some onset temperature, 
T n , where the dynamics starts becoming heterogeneous 
[llGil. The same is true for the (2)-TLG where the SE 
violation is particularly precipitous: D s ~ below T , 
where £ « 0.58 Strong glassformers, on the other 

hand, are anticipated to have only a weak SE violation 
in 2 and 3 dimensions 16] . This is again true for the 
(l)-TLG where, depending on the definition of t q , £ is 
between 0.88 and 1 pi]. 

Therefore, for the interpolated crossover model, as T 
is lowered, DgT a may vary non-monotonically. Indeed, 
this behavior is seen in FIG. [3 Moreover, the smaller 
the value of 6, the larger the extremum. Since 6 varies 
in the same way as T x , this implies that systems with 




FIG. 2: The product of the self-diffusion constant, Ds, and 
structural relaxation time, T a , as a function of temperature 
at various values of the asymmetry parameter, b. Symbols as 
in FIG. □ 



lower crossover temperatures will exhibit a larger decou- 
pling extremum during a fragile-to-strong transition. Of 
course, this trend is only strictly true if strong glassform- 
ers did not themselves violate SE, which is approximately 
true for strong glassformers in 2 and 3 dimensions. The 
analysis of Ref. [6j suggests that salol has a fragile to 
strong crossover at T x « 1.1 Tg. It is at this temper- 
ature where salol also shows non-monotonic decoupling 
|ll|. In the models we study here, we have the same coin- 
cidence of fragile-to-strong crossover and non-monotonic 
decoupling. 

Finally, we present results on the growth of dy- 
namic heterogeneity lengthscales during a fragile-to- 
strong transition. We can study dynamical lengthscales 
by observing a trajectory over a time At and consid- 
ering the binary field n r (t;At) = p r (t)[l — p r (t + At)]. 
Here, p T (t) is 1 if there is a particle at lattice site r 
at time t and otherwise. The field, n T (t;At), gives 
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FIG. 3: Growth of dynamical heterogeneity length as a func- 
tion of structural relaxation time for different values of b. 
Symbols as in FIG.[T] 
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a signal only when there is particle motion at r over 
the range of time between t and t + At. We then ex- 
tract a dynamical lengthscale, £(At), as the inverse of 
the first moment of the structure factors of the mobile 
particles |g, Il2l l2Cl l2ll |22j . These structure factors are 
four point functions: they measure a correlation func- 
tion which depends on two points in time, t and t + At, 
and two points in space, r and r'. The growth of the 
dynamic heterogeneity lengthscale with structural relax- 
ation time also exhibits crossover behavior. At some 
crossover lengthscale, £ x , the growth changes from fragile 
growth to strong growth. This behavior is consistent with 
|6j , where it is argued that generally the relaxation of all 
fragile systems will obey strong dynamics at large enough 
lengthscales. We predict, at least for systems exhibiting 



entropic crossovers 23] like the one studied here, that 
lower crossover temperatures, T x , will give rise to larger 
decoupling extremum as temperature is lowered through 
the crossover region. 
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